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We study the effects of hadronic rescattering on hadron distributions in high-energy nuclear col¬ 
lisions by using an integrated dynamical approach. This approach is based on a hybrid model 
combining (3-|-l)-dimensional ideal hydrodynamics for the quark gluon plasma (QGP), and a trans¬ 
port model for the hadron resonance gas. Since the hadron distributions are the result of the 
entire expansion history of the system, understanding the QGP properties requires investigating 
how rescattering during the hadronic stage affects the final distributions of hadrons. We include 
multistrange hadrons in our study, and quantify the effects of hadronic rescattering on their mean 
transverse momenta and elliptic flow. We find that multistrange hadrons scatter less during the 
hadronic stage than non-strange particles, and thus their distributions reflect the properties of the 
system in an earlier stage than the distributions of non-strange particles. 

PACS numbers: 25.75.Ld, 12.38.Mh, 24.10.Nz, 25.75.Dw 


I. INTRODUCTION 

Quark gluon plasma (QGP), strongly interacting mat¬ 
ter composed of quarks and gluons, has been created in 
high-energy nuclear collisions at the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Labora¬ 
tory (BNL) and the Large Hadron Collider (LHC) at 
CERN [I[. One of the main discoveries at the RHIC 
was that the QGP behaves like a nearly perfect fluid 
01, and since then the transport properties of the QGP, 
especially its shear and bulk viscosity coefficients, have 
been under active investigation. 

The QGP created in high-energy nuclear collisions ex¬ 
pands, cools down, and turns into a hadron gas. Since 
hadrons scatter with each other via the strong interac¬ 
tion, information about the QGP stage is, in general, 
contaminated by rescatterings in the hadronic stage. To 
probe the QGP more directly, thermal photons and dilep¬ 
tons were proposed as penetrating probes [ 1 , 0 . Once 
emitted, photons and dileptons propagate to detectors 
without rescattering since they interact only electromag- 
netically and their mean free path is therefore longer 
than the typical size of the system. However, photons 
and dileptons are not perfect probes of the QGP either. 
They are emitted during all the stages of high-energy 
nuclear collisions such as the primary hard scatterings, 
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pre-equilibrium stage, and late hadronic stage 0 , and 
it is not possible to distinguish photons and dileptons 
emitted from the QGP from those originating from the 
other processes. 

In this paper, we follow up the previous studies [nt 
[13 about multistrange hadrons, and claim that their 
distributions allow to probe the QGP immediately af¬ 
ter hadronization. Multistrange hadrons, in particular 
(j) mesons and Ul baryons, have small scattering cross- 
sections with pions since they do not form any reso¬ 
nances unlike other hadrons. This is a unique property of 
the multistrange hadrons: Their distributions reflect the 
properties of the system mainly at a specific stage deep 
inside the fireball, unlike photons and dileptons which 
come from the all stages. 

By using a hybrid model in which hydrodynamic de¬ 
scription of the QGP fluid is followed by a hadronic cas¬ 
cade model, violation of the mass ordering of the differ¬ 
ential elliptic flow parameter V 2 {pt) was predicted [l5l| . 
This phenomenon was observed recently by the STAR 
Gollaboration [l3, which indicates that (j) mesons are 
less affected by hadronic rescatterings and, consequently, 
that they are good probes deep inside the matter. In 
this paper, we investigate the hadronic rescattering ef¬ 
fects on observables more systematically by employing 
an integrated dynamical approach which we have previ¬ 
ously used to analyze various observables at the RHIG 
and the LHC energies In the previous study of the 
violation of mass ordering [l3 , a first order phase transi¬ 
tion model for the equation of state (EoS) was employed 
and only hadrons up to the mass of A(1232) were taken 
into account in the hadron phase. In the present paper, 
we use a more realistic equation of state |20| connecting 
a parametrized lattice QCD EoS at large temperatures 
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to a hadron resonance gas EoS at low temperatures. The 
hadron resonance gas in the EoS contains all the reso¬ 
nances in the hadron cascade model JAM (Jet AA Mi¬ 
croscopic transport model) [2lj used as the last stage of 
the integrated dynamical approach. Hence, in addition 
to (/), we are able to investigate hadronic rescattering ef¬ 
fects on heavier multistrange hadrons such as S and H. 

In the following we introduce the integrated dynami¬ 
cal model in Sec. im We compare transverse momentum 
(pt) spectra and pT-differential elliptic flow parameters 
V 2 {pt) with the STAR data in Sec. IIIII After we confirm 
we reasonably reproduce the data by using this model, 
we discuss the hadronic rescattering effect in Sec. IIVI In 
this section, we first revisit the violation of mass order¬ 
ing of V 2 (pt)- We next focus on the mean transverse 
momentum (pt) and pT-averaged V 2 calculated with or 
without hadronic rescatterings and investigate how much 
these final observables reflect the information when the 
hydrodynamic stage finishes. We summarize our results 
in Sec. [V] 


II. THE MODEL 

We describe the space-time evolution of high-energy 
nuclear collisions by an integrated dynamical approach 
on an event-by-event basis. We divide the whole reaction 
into three separate stages: Initial stage, hydrodynamic 
stage, and transport stage. The initial entropy density 
distribution after the collision of energetic heavy nuclei 
is parametrized by using a Monte-Carlo version of the 
Glauber model. The subsequent expansion of the matter 
is described by relativistic ideal hydrodynamics. When 
the system becomes sufficiently dilute due to strong ex¬ 
pansion, we switch the description of the system from 
hydrodynamics to kinetic theory. 

In the following we briefly overview each part of the 
model and the corresponding interfaces. Eor further de¬ 
tails, see Ref. M- 


A. Hydrodynamics and equation of state 

Hydrodynamics is a macroscopic description and an ef¬ 
fective theory of the system’s long wavelength/time be¬ 
havior. The hydrodynamical equations of motion are ob¬ 
tained from the conservation of energy and momentum: 

= 0 , ( 1 ) 

where is energy-momentum tensor. In the ideal fluid 
approximation, the energy momentum tensor can be de¬ 
composed as 


Tf^‘' = eu^^u'' - ( 2 ) 

where e, P, and = diag(-|-l, —1, —1, —1) are en¬ 
ergy density, pressure, flow velocity, and the Minkowski 


metric, respectively. To close the system of partial dif¬ 
ferential equations we need to know the equation of state 
(EoS) of the fluid to obtain pressure P as a function of 
energy density e. Once the EoS is known, the space-time 
evolution of thermodynamic quantities and flow veloc¬ 
ity is determined for given initial conditions, i.e., energy 
density distribution and flow field at initial time tq. 

We ignore the baryon number and its continuity equa¬ 
tion since at the collider energies, net baryon density at 
midrapidity is tiny. When we solve Eq. o, temperature 
is not needed since it is merely a parameter to connect 
energy density and pressure. Nevertheless, we also follow 
the space-time evolution of temperature to decide where 
to switch from fluid to cascade. We employ the s95p- 
vl.I parametrization for the EoS. The s95p parametriza- 
tion ( 2 ^ connects the lattice QCD based EoS [ 2 ^ in 
the high temperature region to a hadron resonance gas 
EoS in the low temperature region. The hadronic part 
of s95p-YlA contains the same hadronic species as the 
JAM hadronic cascade model [2l| described in the next 
subsection. 

We numerically solve Eq. ([I} in the Milne coordinates 
(r, r]s,x, y) which is appropriate for the description of the 
evolution at relativistic energies [ 2 ^. Here r = y/t'^ — z'^ 
is longitudinal proper time, ps = ^ lii[(t + z)/{t — z)] 
is space-time rapidity, and x and y are transverse coor¬ 
dinates perpendicular to the collision axis. In the nu¬ 
merical calculations, we employ the Piecewise Parabolic 
Method [ 2 ^ which is known to be able to describe shock 
waves. Eor details about this numerical algorithm, see 
also Ref. [ 2 ^. 


B. Particlization and hadronic cascade 


We employ the hadronic cascade model JAM [2lj to 
describe the space-time evolution of hadron gas after 
“particlization”, i.e., after switching from fluid to par¬ 
ticles. At the late stage of collisions, the system is too 
dilute to maintain equilibrium. We assume this hap¬ 
pens around temperature, which we call the switching 
temperature Tgw, and switch description from h ydr o dy- 
namics to kinetic theory on T = Tsw isosurface [1^ \2^ 
1^ . Within the current model, Tsw is rather adjustable 
parameter which controls the final particle ratios. We 
choose Tsw = 155 MeV to reproduce the observed pr 
spectra of pions, kaons, and protons -I- anti-protons in 
the low Pt region [l^ at the full RHIC energy. At 
T(x) = Tsw, we calculate the single particle phase space 
distributions fi{p,x) for all hadrons included in the EoS 
by employing the Cooper-Frye prescription [s^ 


f,{p,x)d^xd^p 


gi _ P ■ Aad^p 

{2Tr)^E exp(p • u/Tsw) ± 1 


(3) 


Here gi is a degeneracy of the hadronic species i, p^ = 
{E,p) is four-momentum of a particle, and Acr^ is a nor¬ 
mal vector of the T{x) = Tsw hypersurface. It is well 
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known that at some momenta the Cooper-Frye formula 
([3]) gives negative contribution for space-like hypersur¬ 
face element (or even time-like hypersurface element with 
negative time component). These negative numbers cor¬ 
respond to in-coming particles, which one cannot treat 
in the hadronic cascade model. So we just neglect these 
contributions and consider only out-going particles. We 
sample particles by going through all the hypersurface 
elements and generate ensembles of hadrons on an event- 
by-event basis. Note that we do not oversample particles 
to gain statistics for each event. This means each hydro- 
dynamic event corresponds to one ensemble of particles. 
Among all the calculations in the integrated dynamical 
approach, sampling particles from particlization hyper¬ 
surface is numerically most expensive, which is crucial in 
event-by-event simulations. In Ref. [l9|, we discussed in 
detail how to sample particles and integrate the Cooper- 
Frye formula ([3]) over momentum at less numerical cost. 

In hadronic cascade models, experimental hadronic 
cross section data are implemented when available. How¬ 
ever, there are many hadronic scattering processes where 
data do not exist. In such a case, we use the additive 
quark model [ssl - l^ to describe the scattering cross sec- 


tion 
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Here is the total nucleon-nucleon cross section, Ui 
is the number of constituent quarks in a hadron, and 
Usi is the number of strange quarks in a hadron. Thus 
hadrons containing strange quark(s) have relatively small 
total cross section in this model. Moreover, <f) mesons and 
n baryons do not form resonances when scattering with 
pions and therefore have much smaller cross sections than 
non-strange hadrons. 

In the simulations, dynamics of the system is described 
by a sum of incoherent two-body collisions. Two particles 
collide with each other if their minimum distance b in the 
center of mass frame is smaller than the distance given 
by the total cross section 

This is a classical, geometrical interpretation of the two- 
body collision cross section. 

Note that in JAM the properties of hadrons are those 
of free particles. The medium modifications of hadron 
properties such as a broadening and mass shift 
have been extensively studied by various further develop¬ 
ments of transport models such as the relativistic Vlasov- 
Uehling-Uhlenbeck (RVUU) model [s^, the hadron- 
string-dynamics (HSD HI) and parton-hadron-string- 
dynamics (PHSD 0 ) transport models, the Giessen 
Boltzmann-Uehling-Uhlenbeck (GiBUU) [4l|, and the 
Isospin Quantum Molecular Dynamics (IQMD) (d^ . 
Nevertheless, since the main difference in the cross sec¬ 
tions of multistrange particles and the other hadrons is 


due to the (non-)existence of resonance channels, we do 
not expect the in-medium modifications to significantly 
affect this difference, and the use of free particle proper¬ 
ties to provide a conservative baseline is justified. 

As mentioned, the interactions in JAM are described 
as two particle scatterings, and multi-particle scatterings 
are not included. It has been argued that multi-pion 
fusion processes in particular would have a si gni ficant 
effect on the yields of proton-antiproton pairs |43|, but 
subsequent calculations have shown that regeneration is 
a ^ 10% level effect at the LHC [3,0. Nevertheless, 
the multi-meson fusion effects for H baryons are much 
smaller than for protons or A [3|j and thus their inclu¬ 
sion would make the later-discussed differences between 
rescatterings of non-strange and multistrange particles 
even larger. 

As a default setting strong interactions and decays 
are simulated in JAM. However, we switch off the decay 
channels for (j) mesons in the present study to investigate 
hadronic rescattering effects on (j) mesons in an efficient 
way. Otherwise, we would need to perform mass recon¬ 
struction of (/) mesons from two kaons after numerical 
simulations. Since the lifetime of (j) mesons (~ 47 fm/c) 


47| is somewhat larger than the typical life time of the 
system (^10 fm/c), the daughter kaons from (j) meson 
decays (or lack of them) are not expected to affect the 
bulk evolution. As we will discuss in Sec. mi the yield of 
A measured by STAR contains the feed down from elec¬ 
tromagnetic decays. This process does not happen in 
the default setting in JAM and, therefore, is simulated 
separately when we analyze A spectra. We also have an 
option in JAM to deactivate all the hadronic rescatter¬ 
ings and/or resonance decays to investigate how these 
affect final observables. 


C. Initial conditions 

Once initial conditions for the hydrodynamic evolu¬ 
tion are fixed, the subsequent evolution is determined, 
and we are able to obtain final particle distributions 
which can be compared with experimental data. The 
pre-thermalization stage between the first contact of col¬ 
liding nuclei and the initial time of hydrodynamic evo¬ 
lution may be described by non-equilibrium field theory, 
although how exactly it could be done is one of the open 
issues in the physics of relativistic heavy ion collisions. 
Thus, we do not even try to describe it and simply rely 
on the Monte-Garlo (MC) Glauber 0 | and the modified 
Brodsky-Gunion-Kuhn (BGK) models [3] to give the 
density distribution at tq in the transverse plane and in 
the longitudinal direction, respectively. The MC Glauber 
model enables us to calculate the number density of par¬ 
ticipants and of binary collisions on an event-by-event 
basis. Our assumption is that the initial entropy density 
is proportional to a linear combination of these two den¬ 
sities. Even at very high collision energies matter profile 
in the longitudinal direction depends on space-time ra- 
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pidity, and boost invariance does not hold. Furthermore, 
in non-central collisions the longitudinal density gradi¬ 
ent at mid space-time rapidity may be non-zero due to 
the different thicknesses of the colliding nuclei at each 
transverse position. This can be deduced from (pseudo- 
)rapidity distribution of hadrons in p(d)-A collisions. In 
the modified BGK model smooth initial longitudinal pro¬ 
file is parametrized by taking account of the difference in 
the local thickness of colliding nuclei. It is noted that we 
do not employ the Monte-Carlo Kharzeev-Levin-Nardi 
model [4^ which was used to calculate initial condi¬ 
tions in our previous hydrodynamic studies since, in 
the present study, we focus on the hadronic rescattering 
effects on transverse dynamics, which are not supposed 
to be sensitive to the choice of initial conditions. 

Assuming local thermalization, energy density distri¬ 
bution e(r = TQ,r]s,x, y) and pressure P(r = To,r]s,x, y) 
are obtained from the entropy density distribution 
through the s95p-Y\.l EoS. Initial flow velocity at tq is 
supposed to be Bjorken’s scaling flow [5l|, namely, = 1 
and = u^‘ = 0. Throughout this study, we fix 

To = 0.6 fm/c. We choose initial parameters to reproduce 
the pt distributions of pions, kaons, and protons -I- anti¬ 
protons in Au-|-Au collisions at = 200 GeV mea¬ 

sured by the PHENIX Collaboration . The centrality 
dependence of multiplicity is controlled by a fraction of 
soft (participants) and hard (binary collisions) compo¬ 
nent in the MC Glauber model, whereas the multiplicity 
in central collisions is fixed by the overall normalization 
parameter of the initial entropy distribution. On the 
other hand, as mentioned particle ratios are controlled 
by the switching temperature. The parameters control¬ 
ling the shape of the rapidity distribution of particles are 
kept the same as in our previous study M- Since here 
we concentrate on observables at midrapidity, we do not 
discuss them further. 

In hydrodynamic simulations, one needs to convert 
the position distribution (a sum of delta functions in 
the transverse plane) of collision points in MC Glauber 
model to a hydrodynamic density distribution. In our 
model we obtain the density at a point x by counting the 
number of collision points or participants within radius 
xq = \/CTin/TT around the point, and divide by the inelas¬ 
tic cross section in p+p collisions at ^/snn = 200 GeV, 
CTin. For further details, see Ref. M- 


III. SPECTRA AND FLOW 

In this section, we compare invariant pT spectra and 
differential elliptic flow parameter V 2 (j>t) of identified 
hadrons with experimental data from the STAR Col¬ 
laboration. In the present study, we analyzed 0.6 mil¬ 
lion “minimum bias” events (defined as Apart > 2 in the 
MC Glauber calculations for hydrodynamic initial condi¬ 
tions). Centrality is defined using an event distribution 
of the charged hadron multiplicity in \irj\ < 0.5 as done 
by the STAR Collaboration [^. Note that error bars in 


the theoretical plots denote the statistical errors. 


A. Transverse momentum spectra 

In Figs. mm and [3l we show px distributions of identi¬ 
fied hadrons around midrapidity in Au-|-Au collisions at 
= 200 GeV. For clarity, the results and the exper¬ 
imental data points are scaled up or down by a common 
factor. 

Figure [T] shows the px spectra of (a) positive pions 
(I77I < 0.5) and (b) positive kaons (I??! < 0.1) compared 
with the STAR data {\y\ < 0.5 for pions and \y\ < 0.1 
for kaons [H). The Jacobian between y and y is taken 
into account to obtain the invariant px spectra in our 
results. Even if we tuned the initial conditions and the 
switching temperature in the model using pion, kaon, and 
proton -I- anti-proton px spectra from the PHENIX Col¬ 
laboration in our previous study, the results are overall 
in good agreement with the STAR data too. 

In particular, the model reproduces the pion data up 
to Pt ~ 3 GeV/c in central (0-12%) collisions. Above 
this, our result gradually deviates from the data due to 
appearance of (semi-)hard components such as recombi¬ 
nation and jet fragmentation. The more peripheral the 
collision, the lower the px where our result begins to de¬ 
viate from the data. Due to the limited px range of the 
data, the same behavior is not seen in the kaon spectra. 
If one looks at the pion spectrum in peripheral collisions 
(60-80% centrality) carefully, one sees that experimental 
data are systematically larger than the results. This is 
because in peripheral collisions there is a small discrep¬ 
ancy between PHENIX and STAR pion yields. 
Nevertheless we did not further tune the parameters to 
reproduce the STAR data because we want to keep the 
same framework as in our previous study la¬ 
in Fig. [21 we compare the px spectra in \r]\ < 0.5 for (a) 
protons and (b) (p mesons with the STAR data [5J, . 

Due to limited statistics, we do not show some data 
points of the proton spectra in high px regions (above 
Pt 4 GeV/c) in 20-40, 40-60, and 60-80% centrality. 
The proton spectra from STAR are corrected for the A 
and E+ feed down. Weak decays do not occur in the 
default setting in JAM, so we are able to compare our 
results directly with the STAR corrected data. Since we 
neglect baryon chemical potential in our model, our re¬ 
sults for protons are slightly below the experimental data 
at some px range. Nevertheless, overall slopes are in good 
agreement with the data below px ^ S GeV/c. On the 
other hand, similar to pions, deviation between the ex¬ 
perimental data and the results gradually increases with 
Pt, and the more peripheral the collision, the lower the 
Pt where this deviation appears. As mentioned in the 
previous section, we switch off the decays of (p mesons 
in the hadronic cascade to be able to analyze (p meson 
spectra directly without resorting to mass reconstruction 
from kaons. Tendencies seen in (p meson spectra are sim¬ 
ilar to those in proton spectra. 
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FIG. 1: (Color Online) Transverse momentum distributions of (a) tt^ (|?;| < 0.5) and (b) (|r;| < 0.1) obtained from the 

integrated dynamical approach (open square) compared with data from the STAR Collaboration [^, (filled square) for 
^/sNN ~ 200 GeV Au+Au collisions. From top to bottom, each spectrum shows the result of 0-12, 10-20, 20-40, 40-60, and 
60-80% centrality multiplied by 10" with n = 0 to —4 for pious, and 0-5, 5-10, 10-20, 20-30, 30-40, 40-50, 50-60, 60-70, and 
70-80% centrality multiplied by 10" with n = 5 to —3 for kaons. 


Figure [3] shows the pr spectra of strange baryons, (a) 
A -I- A {\r]\ < 1.0), (b) S -h S (I 77 I < 0.75), and (c) 
n + n {\r]\ < 0.75), compared with the STAR data [s^ . 
Similar to the proton case, the more peripheral the col¬ 
lisions, the worse the agreement between experimental 
data and the results, especially above px = 1 .5-2.0 
GeV/c. As expected the measured A and S yields are 
larger than the calculated yields. We had chosen the 
switching temperature Tgw to reproduce the observed 
pion to kaon and pion to proton ratios, which leads to a 
lower temperature than the statistical model fit to all ob¬ 
served hadrons [H, , and thus to lower yields of hyper- 

ons. On the other hand it has been argued that the sta¬ 
tistical model with the common chemical freezeout 
temperature does not reproduce the proton and A yields 
simultaneously at the LHC and RHIC energies [^. It 
was claimed |44| that earlier switching from hydrodynam¬ 
ics to UrQMD at energy density ~ 840 MeV/fm^, which 
corresponds to higher temperature than Tgw = 155 MeV, 
would be a key to resolve the issue, but our results do not 
improve signihcantly if we use higher switching temper¬ 
ature Tsw = 165 MeV. It remains to be seen whether the 
origin of this discrepancy lies in a different description of 
the proton-antiproton annihilations in JAM and UrQMD 
cascades, or in different expansion dynamics in collisions 
at RHIC and LHC energies. 

Note that the experimental yield of A in STAR con¬ 


tains the feed down from electromagnetic 11° decays 
(I]° —>• A -I- 7 ) [ 13 . Since this is the dominant mode 
of S]° decay with the branching ratio ~ 100% and cr = 
2.22 X 10“^^ m [13, the feed down correction of it is 
experimentally challenging. We simulate this decay in 
JAM and include this contribution in the final A spectra 
in Fig. [3] (a). Since the mass of is close to that of A, 
the yield of the primordial 11° is expected to be of the 
same order as the yield of the primordial A. The ratio of 
to A at switching temperature can be estimated as 



After correcting for the effects of resonance decay and 
rescatterings, we hnd that the final particle ratio of 
to A is around 0.3. We have included both sources of A 
in our results, but in general this should be kept in mind 
when comparing theoretical results with the data. 

Although the number of data points from STAR is lim¬ 
ited for n baryons, their yields and slopes are consistent 
with our results. It is worth noticing here that the recent 
lattice QCD calculations suggest the existence of reso- 
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FIG. 2: (Color Online) Transverse momentum distributions of (a) protons and (b) (fi mesons in I?;! < 0.5 obtained from the 
integrated dynamical approach (open square) compared with data from the STAR Collaboration (^ . [s^ (filled square) for 
^sjvjv = 200 GeV Au+Au collisions. From top to bottom, each spectrum shows the result of 0-12, 10-20, 20-40, 40-60, and 
60-80% centrality multiplied by 10" with n = 0 to —4 for protons, and 0-10, 10-20, 20-30, 30-40, 40-50, 50-60, 60-70, and 
70-80% centrality multiplied by 10" with n = 4 to —3 for <j) mesons. 


nances in gen eral or strange baryon resonances in 

particular |22| . which have not been discovered yet. 

Such resonances, whether they are of Hagedorn 
kind [i^, or those predicted by quark models (see 
e.g. Refs. [H [^), would contribute to the yields of 
strange hadrons, and might somewhat change the dy¬ 
namics of the late hadronic stage (b^ - fTOl l . Some of those 
resonances could be formed in the scattering of </> mesons 
or baryons, but it is unknown whether such resonances 
would significantly affect the average scattering cross sec¬ 
tion of (f) mesons and 17 baryons in temperatures below 
the switching temperature. Investigation of this possibil¬ 
ity is beyond the scope of the present paper. 


B. Elliptic flow 

In Fig. |4] the two-particle cumulant V 2 {pt) of pions, 
kaons, and protons in the midrapidity region (|? 7 | < 1.0) 
are compared with STAR U2{2} data at three cen¬ 
tralities (10-20%, 30-40%, and 50-60%). For pions, we 


^ It is, however, unknown whether the difference between the lat¬ 
tice QCD and hadron resonance gas model trace anomalies indi¬ 
cates the inadequacy of the hadron resonance gas model l64l , or 
the existence of so far undiscovered resonance states [^ . 


reasonably reproduce the experimental data for each cen¬ 
trality class. Due to the limited number of events, it is 
hard to obtain U 2 { 2 } with smaller statistical errors, and 
especially difficult this would be for protons. Thus some 
points with large error bars, in particular in low px re¬ 
gion, are not shown in these figures for clarity. We also 
calculated z; 2 {RP} (not shown), namely V 2 with respect 
to theoretically-known reaction plane, to reduce the sta¬ 
tistical errors and found that U 2 {RP} of pions, kaons, 
and protons from our model are in reasonable agreement 
with the STAR data in mid-central collisions. However, 
U 2 {RP} lacks initial fluctuations of the event plane angle 
and, consequently, is significantly smaller than the data 
in both central and peripheral collisions. It is noted that, 
from the MC Glauber model analysis, the fluctuation ef¬ 
fect on eccentricity is large in central and peripheral colli¬ 
sions, but almost negligible in mid-central collisions. See, 
e.g., Fig. 3 (a) in Ref. [t^. At the centralities shown in 
Fig.|3]a clear mass ordering behavior among pions, kaons, 
and protons, namely larger V 2 for smaller mass at fixed 
Pt, is seen both in our results and in the STAR data. 

In Fig. [5l we show the event plane V 2 {pt) of (j), A, 
S, and D at some centrality classes from the integrated 
dynamical approach. In comparison, STAR data from 
the event plane method U 2 {EP} are also shown 

in the figure. We calculate U 2 {EP} by employing the p- 
subevent method using hadrons in \p\ < 1.0 as done by 
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FIG. 3: (Color Online) Transverse momentum distributions of (a) A (|r;| < 1.0), (b) H (|t;| < 0.75), and (c) fl (|n| < 0.75) 
obtained from the integrated dynamical approach (open square) compared with data from the STAR Collaboration (filled 
square) for ^sjvat = 200 GeV Au-fAu collisions. From top to bottom, each spectrum shows data of 0-5, 10-20, 20-40, 40-60, 
and 60-80% centrality multiplied by 10" with n = 0 to —4 for A and S and 0-5, 20-40, and 40-60% centrality multiplied by 10" 
with n — 0 to —2 for fi. Note that the pr distributions of A include the contribution from decay, see the text for details. 
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FIG. 4: (Color Online) Transverse momentum dependence of elliptic flow parameter n2{2} of pions (open sqnare), kaons (open 
circle), and protons (open triangle) in |r;| < 1.0 obtained from the integrated dynamical approach compared with V 2 { 2 } data 
from the STAR Collaboration (filled symbol) for ^/snn = 200 GeV Au+Au collisions. Centrality classes are (a) 10-20%, 
(b) 30-40%, and (c) 50-60%. 


the STAR Collaboration (t^ . 


t;2{EP} = ^(cos[2(<^± - vI/2.^)]). (7) 


The event plane angle '^ 2 ,+ ('I' 2 ,-) is defined for the 
particles with positive (negative) pseudorapidity. In the 
actual experimental analysis, an pseudorapidity gap of 
IA 77 I = 0.075 between two subevents was introduced to 
suppress the nonflow effects [t^. We found, however, 
that introduction of the gap does not change final results 
within error bars in our analysis. The resolution factor 
in this case is calculated as 77. = (cos[ 2 (T' 2 ,+ — 'I' 2 ,-)]). 

In calculating the Q vector, 


Q = ^Wj exp(2i(/)j), ( 8 ) 

3 


needed for the event plane angle 


^' 2 ,± = 



(9) 


a weight factor Wj = mhi{pT,j, 2 GeV/c} was also intro¬ 
duced [t^. However, within statistical errors the results 
are similar to those obtained using Wj = 1. Therefore, 
we employ Wj = 1 in all the results shown in Fig. [5l 
For A and S our results are systematically larger than 
the STAR data. For (p and H our results show similar 
behavior, but since the data have large errors, no firm 
conclusion can be made. 

Our results of elliptic flow parameters for identified 
hadrons are systematically larger than the STAR data. 


which suggests that there is room for finite, but perhaps 
small, viscosity in the fluid-dynamical stage. However, 
we do not go into details about the QGP viscosity since 
the purpose of the present study is to investigate the 
effects of rescattering in the late hadronic stage, which 
are not affected by the QGP viscosity nor by the initial 
state of the system. 


IV. HADRONIC RESCATTERING EFFECTS 


In this section, we investigate hadronic rescattering 
effects on final observables, in particular, how much fi¬ 
nal observables reflect the properties of the system at 
particlization in the integrated dynamical approach. If 
some observables do not change much during the late 
kinetic stage, these observables can be utilized as “pen¬ 
etrating” probes. Namely, their distributions reflect the 
properties of the system immediately after hadroniza- 
tion, and are not contaminated by later hadronic rescat¬ 
terings. To quantify the effects of hadronic rescatter¬ 
ing we perform simulations of nuclear collisions with the 
following three options in the hadronic cascade calcula¬ 
tions: (I) full dynamical evolution (default), (H) deacti¬ 
vating hadronic rescatterings, and (HI) deactivating both 
hadronic rescatterings and resonance decays. To simplify 
the calculations we calculate elliptic flow parameters with 
respect to the theoretically-known reaction plane. 

In Figs, m (a) and (b), we show V 2 [pt) of pions, kaons, 
protons, and (p mesons in minimum bias Au-|-Au colli- 
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FIG. 5: (Color Online) Transverse momentum dependence of the elliptic flow parameter t; 2 {EP} of (a) 0 mesons (0-5, 10-40, 
and 40-80% centrality), (b) A baryons (0-10, 10-40, and 40-80% centrality), (c) H baryons (0-10, 10-40, and 40-80% centrality), 
and (d) H baryons (0-10 and 0-80% centrality) in Ir?! < 1.0 obtained from the integrated dynamical approach (open symbol) 
compared with data from the STAR Collaboration [H^. It^I (filled symbol) for y/SNN = 200 GeV Au+Au collisions. 


sions at ,/snn = 200 GeV from the integrated dynami¬ 
cal approach. Figure IH (a) represents the results without 
hadronic rescatterings (option (II)). These results exhibit 
the mass ordering behavior, namely t;J(pT) > 'c^(pt) > 
v^iPr) > V 2 {pt) for < rriK < rrip < m^, in the 
low pt region. In general the mass ordering results from 
collective flow in which all the components of the fluid 
flow at a common fluid velocity [t^. On the other hand, 


Fig. ini (b) shows a violation of mass ordering between 
protons and </> mesons below pT ~ 1.5 GeV/c: V 2 < V 2 
even though Wp < m^. If one compares these two plots, 
V 2 {pt) decreases below ^1.5 GeV/c during the hadronic 
rescattering stage. On the other hand, v^ipr) stays al¬ 
most unchanged during the hadronic rescattering stage. 
As a consequence, the order is inversed between protons 
and (f) mesons. To see this behavior more clearly, we also 
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show the ratio v^jv^. Below 1 GeV/c, this ratio becomes 
clearly positive indicating the violation of mass ordering. 

The violation of mass ordering was predicted in 
Ref. Ila a nd observed recently by the STAR Collabo¬ 
ration [1^. It should be noted that a larger switching 
temperature (Tsw = 169 MeV) and a first order phase 
transition in the EoS were employed in the calculations 
of Ref. [l^. We confirm in the present study that the 
violation of mass ordering is also observed when a more 
realistic EoS from lattice QCD and a lower switching 
temperature (Tg^ = 155 MeV) are used. 

We can interpret the violation of mass ordering as a 
result of evolution of mean pT and py-averaged V 2 during 
the late transport stage. The slope of the pT-differential 
V 2 can be approximated by the ratio of pr-averaged V 2 
to mean px [7| , 

dV2{pT) _ V2 , . 

dpT (pt) ■ ^ ^ 


If V 2 {pt) is linearly proportional to pr, it is easy to show 
that Eq. (ITUl) holds exactly. Thus a change of slope in 
V 2 {pt) due to hadronic rescatterings could be attributed 
to changes of V2 and/or (pt), which provides us with 
more intuitive understanding. 

We next calculate (px) and p-r-averaged V 2 for each 
identified hadron in |p| < 1 in minimum bias Au-|-Au 
collisions at ^/snn = 200 GeV with the three different 
options mentioned above. We show (px) and V 2 as func¬ 
tions of hadron mass in Figs. [7] (a) and (b), respectively. 
Default results using the option (I) are shown with closed 
squares. To see the effect of hadronic rescatterings, we 
use the option (II) in the hadronic cascade calculation. 
We also use the option (III) to obtain the results ex¬ 
cluding the contributions of both hadronic rescatterings 
and resonance decays. In the default option (I) in which 
both hadronic rescatterings and resonance decays are in¬ 
cluded, (px) is not a linear function of the hadron mass. 
By switching off hadronic rescatterings by employing the 
option (II), we are able to extract the hadronic rescatter¬ 
ing effects on evolution of mean px in the late hadronic 
stage. The difference of (px) between with and without 
hadronic rescatterings are relatively smaller for tt, </, S, 
and n than for the others. In order to confirm a pure flow 
effect we also calculate (px) with option (III) in which 
resonance decay contributions are also not included. As 
is well-known in analysis using the blast wave model [t^ , 
(px ) is a linear function of the hadron mass in this case 
[T^ . We also confirm that it is the case even in our fully 
dynamical calculations. pT-averaged V 2 for three options 
(I), (II), and (III) are compared with each other in Fig. [7] 
(b). By comparing the results (I) and (II) one finds that 
due to hadronic rescatterings V 2 for hadrons other than 
pions increases little - by at most ~6% - but V 2 for pions 
increases by ~25 even though elliptic flow parameters are 
commonly believed to be sensitive to the early stage of 
the reaction. On the other hand, the effect of hadronic 
rescattering on pion mean px is almost absent Q ■ It is 
interesting to note that the effect of resonance decays on 


pion V 2 is also significant, see Fig. [7](b). How resonance 
decays reduce pion V 2 was discussed in Ref. As we 
expected, neither (px) nor V 2 for multistrange hadrons 
is affected by hadronic rescatterings, but at least one of 
these observables is affected for all other particles. 

To see how much hadronic rescatterings affect the final 
observables, we show ratios of (px) and V 2 at particliza- 
tion but after decays (option (II)) to the ones in the final 
state (option (I)) for each identified hadron in Fig.|8l 
For comparison, we also show the ratio obtained from a 
mx scaling ansatz (open square) [t^ in the case of mean 
Px- In this ansatz, one can parametrize px distribution 
as 


dN 

pxdpx 


TeS = 


mx = 


exp 


mx 

Feff 


Tf + -mvf, 


Px 


m^ 


( 11 ) 

( 12 ) 

(13) 


where Tj, m^ and Vf are freeze-out temperature, particle 
mass, and flow velocity, respectively. Parameters were 
chosen to reproduce the slope of the proton px spectrum 
in option (II) (Tf = 0.155 GeV and Vf = 0.56) and in 
option (I) (Tf = 0.130 GeV and Vf = 0.71). If all the 
particles interact strongly with each other and move at 
a common collective velocity during the hadronic stage, 
mass ordering would appear. In fact, the result from the 
mx scaling depends monotonically on particle mass, and 
pions, kaons, and protons follow this monotonic tendency 
obtained from mx scaling ansatz. However, multistrange 
hadrons (including “hidden” strangeness) (</>, S, and H) 
in the integrated dynamical approach deviate from the 
pattern. This means that multistrange hadrons freeze 
out earlier than other hadrons and do not participate in 
the radial flow during the transport stage. Multistrange 
hadrons couple weakly with the system consisting mostly 
of pions since they have small scattering cross sections 
and hardly rescatter with pions in the late hadronic stage. 

Mean transverse momentum (px) for dominant con¬ 
stituents of the medium, namely pions, decreases during 
the hadronic rescattering stage since pdV work done in 
the longitudinal direction reduces the transverse energy 
per unit rapidity (7^.[80l| but the number of pions is fixed. 
In this way, 

(Pt) {Rnid) ^ 0-/ d^Tr) {dEx / dy) ^ ^ (14) 

(T^t) (final) (f/A^ tt) (dF/7^/dy)(final) 

In fact, the two parameter sets in the mx scaling ansatz 
mentioned above are chosen to obey this inequality (fTTl) . 
As for V 2 in Fig. [8] (b), about 20% of final V 2 of pions 
is generated during the late transport stage. In other 
words, only ^80% of final V 2 reflects elliptic flow gen¬ 
erated during the fluid-dynamical stage. Obviously, this 
could depend on centrality (and perhaps collision energy) 
and the number can be regarded as an average value at 
the RHIC energy. Whereas hadronic rescatterings af¬ 
fect little final V 2 for the other hadrons. By combining 
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FIG. 6: (Color Online) Transverse momentum dependence of elliptic flow parameter V 2 for pions (open sqnare), kaons (filled 
square), protons (open circle), and (j) mesons (filled circle) in It;! < 1.0 obtained from the integrated dynamical approach (a) 
without hadronic rescattering and (b) with hadronic rescattering in minimum bias ^sjvjv = 200 GeV Au+Au collisions. The 
lower panels of the plots show the ratio of Wj to uf- 




FIG. 7: (Color Online) (a) Mean transverse momentum (px) and (b) pr-averaged V 2 in minimum bias ^/snn = 200GeV Au+Au 
collisions around midrapidity (|i/| < 1.0) obtained by the hybrid model as a function of the hadron mass. Each point represents 
the result for tt, K, p, (j), A, H, and Q, from left to right. Three lines show the results with different options respectively: (I) 
full calculations with both resonance decays and hadronic rescsatterings (filled square), (II) calculations with resonance decays 
only (open square), and (III) calculations without resonance decays or hadronic rescatterings (open circle). Lines are to guide 
the eye. 
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FIG. 8: (Color Online) Ratio of observables at the fluid stage to the final ones in minimum bias ^^snn = 200 GeV Au+Au 
collisions in \y\ < 1.0 obtained by the integrated dynamical approach (filled square) for (a) mean transverse momentum (pr) 
and (b) pr-averaged V 2 - In comparison, ratio of (pr) obtained from a transverse mass scaling ansatz (open square) is also 
shown (see the text). Each point corresponds to a hadron as shown in the caption of Fig. [7] 


these two results in Fig.|8l we conclude that multistrange 
hadrons are less affected by hadronic rescatterings and, 
therefore, can be used to probe the hadronization stage 
in high-energy nuclear collisions. 

Coming back to violation of mass ordering in differen¬ 
tial V 2 shown in Fig. [6] and also in Fig. 5 in Ref. [H, 
one can interpret this intriguing phenomenon as follows. 
For pions pp-averaged V 2 increases but (pp) changes little 
in the hadronic stage, so slope of uj(pp) gets steeper as 
seen from Eq. m- On the other hand, for protons, pp- 
averaged V 2 does not change so much, but (pp) increases 
in the hadronic stage, so v^ipr) decreases. However, for 
(j) mesons, both pp-averaged V 2 and (pp) after the hydro- 
dynamic stage are almost the same in the final state and, 
consequently, v^ipx) is unchanged. 

Finally we show the normalized freeze-out time (de¬ 
fined as the longitudinal proper time r of the last in¬ 
teraction in hadronic cascade simulations) distributions 
around mid-rapidity (|p| < 1.0) in minimum bias Au-t-Au 
collisions at ^/snn = 200 GeV. In Fig. [9] the results are 
shown in two panels, (a) for mesons and (b) for baryons, 
for clarity. The shaded areas represent the particliza- 
tion time distribution for charged hadrons. Since the 
particlization time does not depend on scattering cross 
sections, it is almost equal for all hadrons within error 
bars. A prominent peak with small width is seen below 
T = 10 fm/c, corresponding to the time of particliza¬ 
tion at different parts of the system. The freeze-out time 
distributions of (p and H do not differ from the particliza¬ 


tion time distribution significantly, which reflects their 
small scattering cross sections. On the other hand, other 
hadrons have broader freeze-out time distributions. For 
S, a peak also exists at an early time like for p and H 
but its height is lower and the peak is accompanied by 
a broad tail. This is because S has a contribution from 
long-living resonance S(1530) (decay width F = 9.9 MeV 
~ 1/(20 fm/c) [13) whose decay forms a long tail in the 
freeze-out time distribution. So the primordial S freezes 
out as early as p and H do. These results support our 
findings above that multistrange hadrons freeze out just 
after the particlization due to less rescatterings in the 
late kinetic stage. 


V. SUMMARY 

We have studied the effects of the hadronic rescat¬ 
terings especially on multistrange hadrons and claimed 
they can be used to probe the hadronization stage in 
high-energy nuclear collisions. We have employed an 
integrated dynamical model in which the ideal hydro- 
dynamic model is combined with the hadronic cascade 
model, JAM, and confirmed that our model fairly repro¬ 
duces the experimental data of pp spectra and differential 
V 2 not only for pions, kaons, and protons but also for mul¬ 
tistrange hadrons. Within this approach, we have simu¬ 
lated the hadronic stage with or without the rescatterings 
and compared the calculated mean transverse momen- 
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FIG. 9: (Color Online) Normalized freeze-out time (r) distributions for (a) mesons: tt, K, and (j), and (b) baryons: p, A, H, and 
Q in \y\ < 1.0 in minimum bias ^snn = 200 GeV Au-|-Au collisions. The shaded areas correspond to normalized particlization 
time distributions of charged hadrons. 


turn and py-averaged V 2 in these two cases to investigate 
the hadronic rescattering effects. We have found that 
the multistrange hadrons are less affected by the rescat¬ 
terings than non-strange hadrons. Because of the small 
scattering cross sections, the multistrange hadrons do not 
fully participate in the radial flow during the hadronic 
stage and freeze out earlier than non-strange hadrons. 
With these results we have shown how to interpret behav¬ 
iors of mass ordering and its violation in V 2 {pt)- Changes 
of slope of V 2 {pt) during the hadronic rescattering stage 
result from interplay between changes of mean pp and 
PT-averaged V 2 ■ We have also showed the freeze-out time 
distributions for identified hadrons and confirmed that 
the multistrange hadrons freeze out soon after they are 
“particlized”. 

In future, we also plan to investigate the hadronic 


rescattering effects on observables for multistrange 
hadrons at the LHC energy. Although elliptic flow pa¬ 
rameters are expected to be less affected by the hadronic 
rescatterings at the LHC, it would be important to quan¬ 
tify how much they do change, and thus how much the 
observed V 2 at LHC reflects the QGP evolution. 


Acknowledgments 

We acknowledge a fruitful discussion with N. Xu. This 
work was supported by JSPS KAKENHI Grant Numbers 
12J08554(K.M.) and 25400269 (T.H.), and by BMBF un¬ 
der contract no. 06FY9092(P.H.). 


[1] K. Yagi, T. Hatsuda, and Y. Miake, Gamb. Monogr. Part. 
Phys. Nucl. Phys. Gosmol. 23, 1 (2005). 

[2] U.W. Heinz and P.F. Kolb, Nucl. Phys. A 702, 269 

( 2002 ). 

[3] M. Gyulassy, arXiv:nucl-th/0403032 

[4] T.D. Lee, Nucl. Phys. A 750, 1 (2005). 

[5] M. Gyulassy and L. McLerran, Nucl. Phys. A 750, 30 
(2005). 

[6] E.V. Shuryak, Nucl. Phys. A 750, 64 (2005). 

[7] T. Hirano and M. Gyulassy, Nucl. Phys. A 769, 71 
(2006). 

[8] E.L. Feinberg, Nuovo Cim. A 34, 391 (1976). 


[9] E.V. Shuryak, Phys. Lett. B 78, 150 (1978). 

[10] E.L. Bratkovskaya, Nucl. Phys. A 931, 194 (2014). 

[11] A. Shor, Phys. Rev. Lett. 54, 1122 (1985). 

[12] H. van Hecke, H. Sorge, and N. Xu, Phys. Rev. Lett. 81, 
5764 (1998). 

[13] S.A. Bass and A. Dumitru, Phys. Rev. C 61, 064909 

( 2000 ). 

[14] Y. Cheng, F. Liu, Z. Liu, K. Schweda, and N. Xu, Phys. 
Rev. C 68, 034910 (2003). 

[15] T. Hirano, U. Heinz, D. Kharzeev, R. Lacey, and Y. Nara, 
Phys. Rev. C 77, 044909 (2008). 

[16] M. He, R.J. Fries, and R. Rapp, Phys. Rev. C 85, 044911 





















14 


( 2012 ). 

[17] X. Zhu, F. Meng, H. Song, and Y.X. Liu, Phys. Rev. C 
91, 034904 (2015^ 

[18] M. Nasim [STAR Collaboration], Nucl. Phys. A 904-905, 
413c (2013). 

[19] T. Hirano, P. Huovinen, K. Murase, and Y. Nara, Prog. 
Part. Nucl. Phys. 70, 108 (2013). 

[20] P. Huovinen and P. Petreczky, Nucl. Phys. A 837, 26 

( 2010 ). 

[21] Y. Nara, N. Otuka, A. Ohnishi, K. Niita, and S. Chiba, 
Phys. Rev. C 61 , 024901 (2000). 

[22] A. Bazavov et al, Phys. Rev. Lett. 113 , 072001 (2014). 

[23] T. Hirano, Phys. Rev. C 65, 011901 (2002). 

[24] P. Colella and P.R. Woodward, J. Comput. Phys. 54, 174 
(1984). 

[25] T. Hirano and Y. Nara, Prog. Theor. Exp. Phys. 2012 , 
01A203 ( 2012 ). 

[26] S.A. Bass, A. Dumitru, M. Bleicher, L. Bravina, 
E. Zabrodin, H. Stoecker, and W. Greiner, Phys. Rev. 
C 60 , 021902 (1999). 

[27] D. Teaney, J. Lauret, and E.V. Shuryak, Phys. Rev. Lett. 
86 , 4783 (2001). 

[28] D. Teaney, J. Lauret, and E.V. Shuryak, 

arXiv:nucl-th/0110037 

[29] T. Hirano, U.W. Heinz, D. Kharzeev, R. Lacey, and 

Y. Nara, Phys. Lett. B 636, 299 (2006). 

[30] C. Nonaka and S.A. Bass, Phys. Rev. C 75, 014902 
(2007). 

[31] H. Petersen, J. Steinheimer, G. Burau, M. Bleicher and 
H. Stocker, Phys. Rev. G 78, 044901 (2008). 

[32] F. Gooper and G. Frye, Phys. Rev. D 10, 186 (1974). 

[33] H. Sorge, A. von Keitz, R. Mattiello, H. Stocker, and 
W. Greiner, Z. Phys. G 47, 629 (1990). 

[34] H. Sorge, L. Winckelmann, H. Stocker, and W. Greiner, 

Z. Phys. C 59, 85 (1993). 

[35] S.A. Bass et al, Prog. Part. Nucl. Phys. 41, 225 (1998). 

[36] M. Bleicher et al., J. Phys. G 25, 1859 (1999). 

[37] C.M. Ko, V. Koch, and G.Q. Li, Ann. Rev. Nucl. Part. 
Sci. 47, 505 (1997). 

[38] X.S. Fang, G.M. Ko, G.Q. Li, and Y.M. Zheng, Nucl. 
Phys. A 575, 766 (1994). 

[39] W. Gassing and E.L. Bratkovskaya, Phys. Rept. 308 , 65 
(1999). 

[40] W. Gassing, Eur. Phys. J. ST 168, 3 (2009). 

[41] O. Buss, T. Gaitanos, K. Gallmeister, H. van Hees, 
M. Kaskulov, O. Lalakulich, A.B. Larionov, T. Leitner, 
J. Weil, and U. Mosel, Phys. Rept. 512, 1 (2012). 

[42] C. Hartnack, H. Oeschler, Y. Leifels, E.L. Bratkovskaya, 
and J. Aichelin, Phys. Rept. 510, 119 (2012). 

[43] R. Rapp and E.V. Shuryak, Phys. Rev. Lett. 86 , 2980 

( 2001 ). 

[44] J. Steinheimer, J. Aichelin, and M. Bleicher, Phys. Rev. 
Lett. 110 , 042501 (2013). 

[45] Y. Pan and S. Pratt, Phys. Rev. C 89, 044911 (2014). 

[46] P. Huovinen and J.I. Kapusta, Phys. Rev. C 69, 014902 
(2004). 

[47] K.A. Olive et al. [Particle Data Group[, Chin. Phys. C, 
38, 090001 (2014). 

[48] M.L. Miller, K. Reygers, S.J. Sanders, and P. Steinberg, 
Ann. Rev. Nucl. Part. Sci. 57, 205 (2007). 

[49] H.-J. Drescher and Y. Nara, Phys. Rev. C 75, 034905 


(2007). 

[50] H.-J. Drescher and Y. Nara, Phys. Rev. C 76, 041903 
(2007). 

[51] J.D. Bjorken, Phys. Rev. D 27, 140 (1983). 

[52] S.S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 
69 , 034909 (2004). 

[53] C. Adler et al. [STAR Collaboration[, Phys. Rev. Lett. 
89, 202301 (2002). 

[54] B.I. Abelev et al. [STAR Collaboration[, Phys. Rev. Lett. 

97 , 152301 (2006). 

[55] J. Adams et al. [STAR Collaboration[, Phys. Rev. Lett. 
92, 112301 (2004). 

[56] B.I. Abelev et al. [STAR Collaboration], Phys. Rev. C 
79 , 064903 (2009). 

[57] J. Adams et al. [STAR Collaboration[, Phys. Rev. Lett. 

98, 062301 (2007). 

[58] J. Cleymans, B. Kampfer, M. Kaneta, S. Wheaton, and 
N. Xu, Phys. Rev. C 71 , 054901 (2005). 

[59] P. Huovinen, Eur. Phys. J. A 37, 121 (2008). 

[60] A. Andronic, P. Braun-Munzinger, and J. Stachel, Phys. 
Lett. B 673 , 142 (2009) [Erratum-ibid. B 678 , 516 
(2009)]. 

[61] L. Milano [ALICE Collaboration], Nucl. Phys. A 904- 
905, 531c (2013). 

[62] C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 
89 , 092301 (2002). 

[63] A. Majumder and B. Muller, Phys. Rev. Lett. 105 , 
252002 (2010). 

[64] A. Wiranata, V. Koch, M. Prakash and X.N. Wang, 
Phys. Rev. C 88, no. 4, 044917 (2013). 

[65] R. Hagedorn, Nuovo Cim. Suppl. 3, 147 (1965). 

[66] U. Loring, B. C. Metsch and H.R. Retry, Eur. Phys. J. A 
10 , 395 (2001). 

[67] U. Loring, B. C. Metsch and H.R. Retry, Eur. Phys. J. A 
10 , 447 (2001). 

[68] C. Greiner, P. Koch-Steinheimer, F.M. Liu, 

LA. Shovkovy and H. Stoecker, J. Phys. G 31 , 
S725 (2005). 

[69] J. Noronha-Hostler, J. Noronha, G.S. Denicol, 
R.P.G. Andrade, F. Grassi and C. Greiner, Phys. 
Rev. C 89, no. 5, 054904 (2014). 

[70] J. Noronha-Hostler and C. Greiner, arXiv: 1405.7298 
[nucl-th]. 

[71] J. Adams et al. [STAR Gollaboration[, Phys. Rev. C 72, 
014904 (2005). 

[72] T. Hirano and Y. Nara, Phys. Rev. C 79, 064904 (2009). 

[73] B.I. Abelev et al. [STAR Collaboration], Phys. Rev. C 
77 , 054901 (2008). 

[74] P. Huovinen, P.F. Kolb, U.W. Heinz, P.V. Ruuskanen, 
and S.A. Voloshin, Phys. Lett. B 503, 58 (2001). 

[75] E. Schnedermann, J. Sollfrank, and U.W. Heinz, Phys. 
Rev. C 48 , 2462 (1993). 

[76] H. Heiselberg and A.M. Levy, Phys. Rev. C 59, 2716 
(1999). 

[77] T. Hirano, Phys. Rev. Lett. 86, 2754 (2001). 

[78] T. Hirano, P. Huovinen, and Y. Nara, Phys. Rev. C 83 , 
021902 (2011). 

[79] M. Gyulassy and T. Matsui, Phys. Rev. D 29, 419 (1984). 

[80] P.V. Ruuskanen, Phys. Lett. B 147 , 465 (1984). 



